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ABSTRACT 

We analyse the stellar and hot gas content of 18 nearby, low-mass galaxy clusters, detected in 
redshift space and selected to have a dynamical mass 3 x 10^'* < M/Mq < 6 x 10^'' {h = 
0.7), as measured from the 2dF Galaxy Redshift Survey. We combine X-ray measurements 
from both Chandra and XMM with ground-based near-infrared observations from CTIO, AAT 
and CFHT to compare the mass in hot gas and stars to the dynamical mass and state of the 
clusters. Only 13 of the clusters are detected in X-ray emission, and for these systems we 
find that a range of 7-20 per cent of their baryonic mass, and < 3 per cent of their dynamical 
mass, is detected in starlight, similar to what is observed in more massive clusters. In contrast, 
the five undetected clusters are underluminous in X-ray emission, by up to a factor 10, given 
their stellar mass. Although the velocity distribution of cluster members in these systems is 
indistinguishable from a Gaussian, all show subtle signs of being unrelaxed: either they lack 
a central, dominant galaxy, or the bright galaxy distribution is less concentrated and/or more 
elongated than the rest of the sample. Thus we conclude that low-mass clusters and groups 
selected from the velocity distribution of their galaxies exhibit a dichotomy in their hot gas 
properties. Either they are detected in X-ray, in which case they generally lie on the usual 
scaling relations, or they are completely undetected in X-ray emission. The non-detections 
may be partly related to the apparently young dynamical state of the clusters, but it remains a 
distinct possibility that some of these systems are exceptionally devoid of hot emitting gas as 
the result of its expulsion or rarefaction. 
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1 INTRODUCTION 

Galaxy clusters and groups are important both as cosmological 
probes, and as laboratories for studying galaxy evolution. In par- 
ticular, their deep gravitational potential means that their baryon 
content should be nearly representative of the Universe as a whole, 
and that the diffuse gas is at a temperature that is accessible to ob- 
servation. As a result, they represent one of the few places where it 
is possible to study the stars, cold and hot gas, and dark matter in a 
single system. 

It is now well known that the mass fraction in stars is not 
univ ersal, but in general decreases with increas i ng cluster mass 
(e.g. 'Eke et al. '2004hl l2005l :l RameUa et d] |2004 iLin et alJllOOsI : 
lOiodin i et al. 2009). On the other h and, the mass fraction of hot gas 
appears to increase with mass (e.g. IVikhlinin et alj2006l : ISun et all 
l2009l : |Prattetal.|[2O10l) . If these systems are closed boxes, then 
the sum of stellar and gas mass fractions should equal the univer- 
sal value fb — Qb/^m, where fib and Sim are the baryon and 



matter densities, respectively , relative to the critical density. Re- 
cently iGonzalez et alj ( l2007h claimed this to be the case, in par- 
ticular arguing that in the lowest mass systems there is signifi- 
cant stellar mass in the intracluster-light and halo of the central 
galaxy, which completes the baryon fraction. The interpretation 
then is that most of the baryons in these low-mass groups have 
cooled to form stars; this poses a challenge for normal hierarchi- 
cal models which predic t that more massive clusters are actually 
built from these groups ( lBaloghetal1l2008l) . An important open 
question then is whether or not there really exists a significant pop- 
ulation of groups with A/stars / A/gas > 20 per cent. The conclu- 
sions of iGonzalez et al.l ( |2007|) depend partly on an extrapolated 
mean relation for the gas fraction of clusters as a function of mass, 
taken from Vikhlinin et al. ( 2006). However, similar conclusions 
were reached by Lagana et al. ( 2008), based on a small sample of 
five nearby Abell clusters observed with XMM. 

The other possible explanation for the high ratio of stellar- 
to-gas mass in groups is that these systems are deficient in X- 
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ray emitting g as, an d ass ociated metals. Recent observ ations by 
iGiodini et alj ( |2009|) and iRasmussen & PonmanI ( l2009h indicate 
this is the case; in particular the latter provides an interesting anal- 

i 'sis of the abundances in a sample of 15 X-ray bright groups 
Rasmussen & Ponmanl2007h and concludes that these groups have 
fewer metals than expected given their stellar mass. However, their 
samp le was selected (from the larger sample of lOsmond & PonmanI 
|2004) to have relaxed X-ray morphologies and good photon statis- 
tics. They may therefore be biased toward undisturbed, cool-core 
groups which are unlikely to be the typical precursor of more mas- 
sive clusters I Rasmussen & Ponmanll200^ . The apparent require- 
ment for strong energy sources to counter the high cooling rates 
near the centre of clusters may provide a natural mechanism for 
removing hot gas from groups. Models in corporating supernova- 
driven superwinds (e.g. Dave et al]|2008ah or supermassive black 
hole accretion (e.g. lBower et al.l2008l ; lMcCarthv et al.l2010ah have 
had considerable success at matching the observed properties of 
X-ray emitting gas in groups and clusters. Alternatively, models 
in which gas was heated prior to the virialization of the group 
or cluster predict that such systems would fail to accrete their 
full c omplement of gas in the first place (e.g. iBalogh et al.ll 19991 : 
[McC arthy et al. 2008). 

Most of the work that has been done so far has been based 
on X-ray selected samples of clusters, which introduces a po- 
tential bias. Several investigations have been conducted in an 
attempt to evaluat e the importan c e of t his bias. In particular, 
Bower et alj (Il994h lOilbank et alj j2()04h, lOonahue e t al. ( 2001), 
Popesso et al. I ( l2007h and lRozo et alj make use of pointed 



and archival ROSAT observations of cluster samples, identified 
based upon their optical galaxy population using redshift and po- 
sition information. All of these independent studies found signif- 
icant scatter between the optical richness (or total optical lumi- 
nosity) and X-ray lumi nosity of their system s. In particular, both 
iBower et all ( 1 1 997b and iRvkoff et all l l2008ij) find their optically- 
detected clusters are systematically underluminous in X-rays, rel- 
ative to X -ray selected samp les. Similar conclusions have been 
reached by iHicks et al.l ( |2008|) . using Chandra observations of 13 
optically-selected clusters at 0.6 < z < 1.1. This is at least 
partly due to a Malmquist-type bias in X-ray flu x limited sur - 
veys jlkebeeta l. 2002: Pratt etal. 2009; Vikhlini netal] |2009|) ; 
for example. iRvkoff et alj ( l2008bl) use weak-lensing based mass 
estimates to conclude that, once a correction is made for non- 
hydrostatic equilibrium and this Malmquist bias, the Lx — M rela- 
tion defined by their optically-selected clusters is similar to that of 
the X-ray selected HIFLUGCS ( Reiprich & Bohringer 2002|) sam- 
ple. 

Recently, IRasmussen et alj i2006h have begun an interest- 
ing study of a redshift-selected, unbiased group sample followed 
up with XMM. Of the first nine analysed, only three were de- 
tected, and two of these a re underluminous given their velocity 
dispersion I Bai et al. 20101) . Based on their velocity distribution, 
IRasmussen et alj ^200^ argue that the underluminous systems are 
in the early stage of collapse; similar interpretation s of such sys- 
tems h ave been made by iPopesso et alj ( 120071) and loietrich et alj 
( l2009h . 

What is clear is that the mass and dynamical state distribution 
of a cluster sample will be sensitive to the way the clusters are de- 
tected. It is important to understand these effects, both to identify 
robust indicators of a cluster's mass and to use them as cosmo- 
logical probes. Perhaps even more interesting is the fact that the 
scatter in key properties such as gas entropy and stellar fraction, at 
fixed total mass, contains important information about the physical 



processes associated with galaxy formation (e.g. iMcCarthv et alj 
l2010bl : lDave et alj|2008bl ). 

One of the most common ways to select groups and clus- 
ters from large redshift surveys is the friends-of-friends linking 
method, which has the advantage that it is easy to implement 
and fairly strai ghtforward to calibrate with numerical simulations 



(e.g. lEke et alj 2004a). Eke et al. (2005) used such a sample from 



the 2dFGRS (CoUess et al. 2001) to demonstrate that most of the 
stars in the Universe are associated with small groups, M ~ 
2 X 10^^ Mp), in reasonable agreement with model predictions 
dBower et al TI2OO6I). It would be of interest to know the fate of 



the hot gas in groups selected in this way: in particular, what is 
the efficiency of converting gas to stars on average, and what is its 
variation between systems? 

We have therefore selected 18 low-mass clusters, with 3 x 
10^'* < M/Mq < 6 X 10^*, from the catalogue of lEkeetalJ 
( l2004j) . All but two of the 18 clusters were followed up with 
pointed Chandra or XMM observations. This represents the first 
time that a complete, optically-selected sample of low-mass clus- 
ters has been followed up with these facilities. Surprisingly, we 
found that despite the narrow selection on dynamical mass, the 
clusters have a wide range of X-ray luminosities, and five of them 
were undetected. The details of the object selection and X-ray ob- 
servations are given in Mazzotta et al. (in prep, hereafter Paper II). 
Here we present an analysis of near-infrared (NIR) data in these 
clusters, to measure the stellar luminosity and mass. The plan of 
the paper is as follows. The origin and reduction of the NIR data is 
described in § |2| The details of the photometry and calculation of 
NIR luminosity are given in §|3] We then calculate the total NIR lu- 
minosity of each cluster and consider various scaling relationships 
with dynamical and X-ray properties, in § |4] We find that most 
of the clusters appear normal when compared with X-ray selected 
samples; we explore possible explanations for the five undetected 
clusters in §|5] finally drawing our conclusions in §|6| 

We use a cosmology with Sim = 0.3, f^A = 0.7 and h = 
iyo/(100km/s/Mpc)= 0.7. All NIR magnitudes are on the Vega 
(2MASS) system. 



2 OBSERVATIONS 



2.1 Cluster selection 



We selected 18 clusters from the catalogue of lEke et alj i2004al) ; 
the adopted centres and redshifts, which in some cases differ from 
those in the original catalogue, are given in Table [T] Our main 
aim was to look for the variation in X-ray emission from a mass- 
selected sample; thus we considered all groups from that catalogue 
with 3 X 10" < M/Mq < 6 X lO". Our defa ult definition of th e 
mass, radius and velocity dispersion come from lEke et al. U2004al) . 
where these quantities are related by 



(1) 



We excluded a few groups from this selection with large uncer- 
tainties on the dynamical mass, due to poor membership or clearly 
non-Gaussian velocity dispersions. All selected groups have red- 
shifts for at least 15 members. 

T he velocity disp ersions are computed using the gapper esti- 
mate of lBeers et alj ( ll990i) . and i?rms is the weighted rms projected 
separation from the cluster centre of all members. The factor of 5 
was chosen to give a mass which is in good agreement with dark 
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iCl 


Name 


RA 


Dec 




NIR source 


Image size 






(J2000) 






(ai'cmin) 


1 


A2734 




9Q 801 77 


0.061 


AAT 


21.7 X 21.7 


2 


A3880 


J JO.i/ / JZO 




0.058 


AAT 


21.7 X 21.6 


3 


A3094 


47 QOSSl 


9fi ssns7 


0.068 


2MASS 


39.0 X 55.5 


4 


A1650 




1 80407 


0.084 


CFHT 


27.4 X 27.4 


5 


RBS317 


'ifi 9Q4 


-29 486 


0.060 


AAT 


21.7 X 21.7 


6 


MS 1306.7-0121 


19/.31UD/ 


-1.61012 


0.086 


CTIO 


31.8 X 31.8 


7 


S0041 


6.38084 


-33.04638 


0.050 


2MASS 


38.8 X 65.5 


8 


A954 


153.43704 


-0.12043 


0.095 


CFHT 


27.3 X 27.4 


9 


A1663 


195.60600 


-2.52938 


0.083 


2MASS 


39.6 X 65.4 


10 


Chan4990 


198.03166 


-0.98243 


0.084 


CFHT 


27.3 X 27.4 


11 


XMM5 


48.80986 


-29.14884 


0.068 


2MASS 


38.6 X 55.5 


12 


Chan4991 


170.80439 


1.08753 


0.074 


CTIO 


31.7 X 31.9 


13 


XMM3 


8.96560 


-27.52812 


0.071 


2MASS 


38.8 X 55.4 


14 


XMM9 


201.76353 


1.31428 


0.081 


CFHT 


27.3 X 27.4 


15 


XMM4 


199.83071 


-0.88001 


0.084 


2MASS 


30.0 X 30.0 


16 




150.93404 


-2.17231 


0.096 


CTIO 


31.8 X 31.9 


17 


XMMIO 


202.61707 


1.35255 


0.082 


2MASS 


47.1 X 52.0 


18 




349.50853 


-28.17449 


0.077 


AAT 


21.6 X 21.7 



Table 1. Basic coordinates of all clusters in our sample, including the source of the infrared imaging used in this paper. 



Id 


cr 


Mdyn 




M200 


-R20() 


M500 


-R500 


^max 




(km/s) 


(10" Mq) 


(Mpc) 


(IO^Mq) 


(Mpc) 




(Mpc) 


(R500) 


1 


590±54 


3.6±0.8 


0.89±0.04 


3.25±0.82 


1.34±0.12 


7.06±2.27 


1.32±0.14 


0.64 


2 


580±58 


4.8±1.0 


1.23±0.04 


3.10±0.69 


1.32±0.13 


3.09±1.17 


1.00±0.12 


0.80 


3 


700±53 


5.0±1.0 


0.88±0.04 


5.37±1.25 


1.57±0.12 


1.49±0.42 


0.78±0.07 


1.4 


4 


600±54 


5.6±1.2 


1.34±0.04 


3.31±0.81 


1.32±0.12 


6.94±0.17 


1.30±0.01 


0.71 


5 


550±64 


3.5±1.0 


1.00±0.04 


2.63±0.91 


1.25±0.15 


1.24±0.10 


0.74±0.02 


0.80 


6 


540±69 


3.8±1.2 


1.12±0.07 


2.40±0.82 


1.18±0.15 


1.40±0.09 


0.76±0.02 


0.81 


7 


580±49 


4.1±0.9 


1.05±0.04 


3.13±0.81 


1.33±0.11 


6.04±3.37 


1.21±0.27 


0.26 


8 


680±72 


4.5±1.1 


0.84±0.06 


4.74±1.26 


1.47±0.16 


1.00±0.14 


0.68±0.03 


0.70 


9 


620±63 


5.0±1.1 


1.12±0.04 


3.65±0.96 


1.36±0.14 


1.68±0.02 


0.81 ±0.00 


0.94 


10 


530±70 


3.1±0.9 


0.95±0.05 


2.28±0.72 


1.16±0.15 


2.12±0.88 


0.87±0.12 


1.1 


11* 


570±59 


3.3±0.8 


0.87±0.05 


2.90±0.85 


1.28±0.13 


< 0.18 


< 0.38 




120 


560±57 


3.4±0.8 


0.93±0.04 


2.73±0.77 


1.24±0.13 


0.91±0.45 


0.66±0.10 


0.66 


13* 


640±59 


3.6±1.1 


0.76±0.08 


4.08±1.27 


1.43±0.13 


< 0.18 


< 0.38 




14 


720±61 


5.5±1.2 


0.91±0.06 


5.74±1.32 


1.59±0.13 


0.55±0.07 


0.56±0.02 


0.80 


15 


700±69 


5.6±1.3 


0.98±0.05 


5.25±1.39 


1.53±0.15 


0.74±0.37 


0.61±0.08 


0.54 


16* 


480±111 


3.2±1.3 


1.19±0.09 


1.66±0.77 


1.04±0.24 


< 0.18 


< 0.38 




17* 


730±97 


4.3±1.7 


0.69±0.09 


5.97±2.59 


1.60±0.21 


< 0.18 


< 0.38 




18* 


710±94 


5.2±2.1 


0.89±0.07 


5.53±2.59 


1.57±0.21 


< 0.18 


< 0.38 





Table 2. Dynamical properties, mass and radii measurements for each cluster used in this paper. The radii -Rrms and -R200, with their associated mass 
estimates M^y^ and A/200 are computed from the galaxy redshifts and positions as described in the text. M500 and -R500 are measured from the X-ray 
images, and rmax is the radius within which X-ray emission is detected, in units of rsoo. The starred entries are undetected in X-ray. and for these systems 
we calculate upper limits on -R500 and M500 by extrapolating the correlation between these quantities and Lbol- Groups 16 and 18 do not have Chandra or 
XMM observations. 



matter haloes identified in numerical simulations, using a friends- 
of-friends algorithm with linking parameter b — 0.2. Therefore the 
mass Mdyn corresponds to the mass on scales larger than -Rrms- 
We recalculate uncertainties on a and Rrms using a jackknife tech- 
nique. Uncertainties on a are typically 10-20 per cent, while those 
on the radius are typically 5-10%. These lead to dynamical mass 
estimates which are uncertain by 20-40%. 

For comparison with the literature it will also be useful to de- 
fine a dynamical est imate of the "virial radius", i?200- Following 
lRamellaetan ( l2004h we define 



10//o(l + z)-3/2 



and 

3 2 

M200 = 7^.R20oo- . (3) 

In contrast with Mdyn, this mass depends only on the measured 
velocity dispersion, and is approximately independent of the radial 
extent of the friends-of-friends group, since in general a is a weak 
function of radius outside the core. 

The original 2dFGRS imposed a bright magnitude limit on the 
spectroscopic selection iColless et al. 2001). As a result, the bright- 
est galaxies in several of our clusters do not have a redshift from 
the original survey. Many of these were foun d in the NASA Ex tra- 
galactic Database and in the 6dFGRS survey ( I Jones et al .120091) . In 
three clusters (8, 10, 11) there are still fairly bright galaxies without 
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redshifts, but they are subdominant and it makes little difference 
whether we include them as cluster members or not (we do). For 
cluster 13, the central very bright galaxy has no redshift available, 
and this contributes substantially to the cluster luminosity. As it is 
centrally located we feel confident that it is a cluster member, and 
thus throughout the paper we assume this is the case. 

2.2 X-ray Observations 

The description of the acquisition, reduction and analysis of X- 
ray data from Chandra and XMM are given in Paper II. Here we 
summarize the salient details. Chandra data have been analyzed 
using CIAO v4.0 and CALDB v3.4.3, and the standard tools. 
XMM data have bee n analy zed using the procedure described in 
iBourdin & Mazzott i tarn . For each cluster we extract the sur- 
face brightness and temperature profiles which are then used to es- 
timate th e gas and total mass pr ofile following the approach pro- 
posed by IVikhlinin et alj ( l2006h . This procedure, which assumes 
hydrostatic equilibrium and spherical symmetry, involves model- 
ing the 3-D density and 3-D temperature profiles and fitting the 
projected quantities to the corresponding data set. In the projec- 
tion we take into account the instrumen t responses and the eff ect of 
the "spectroscopic-like" temperatures jMazzotta et alj|2004h . The 
radius and mass at an overdensity of 500, -Rsoo and M500 respec- 
tively, are calculated directly from the estimated cluster total mass 
profile, and are listed in Table[2] Note that because our clusters are 
nearby, low surface brightness systems, most of the X-ray data do 
not extend to ; in Table|2]we list the maximum radius to which 
X-ray emission is detected, in units of i?5oo. For most clusters, this 
is > 70 per cent; but clusters 7 and 15 require significant extrapo- 
lation to obtain Afsoo. The other estimates of mass and radius are 
included in this table, as well. 

Five of the clusters are undetected in X-rays and, for these, 
we estimate an upper limit to the bolometric X-ray luminosity 
of about I/x,boi < 3 X 10*^ergs s^^; although the precise limit 
depends on the unknown system temperature and metallicity, this 
limit is appropriate if kT > 0.2keV (Paper II). This is a con- 
servative limit, corresponding to 1000 counts within the field of 
view. If the clusters were twice as large as this field, we would 
still expect ~ 250 counts; this would be measurable, even in the 
ROSAT data which is all that exists for groups 16 and 18. We will 
therefore adopt this upper limit throughout the paper. Correspond- 
ing upper limits on -R500, Afsoo and Mgas are derived by fitting 
a logarithmic relation between these quantities and Lboi, for the 
13 detected systems, and extrapolating to the upper limit on lumi- 
nosity. This yields i?500 < 0.4 Mpc, M500 < 1.8 x lO^M©, 
and A/gas < 1.3 x 10^^ A/0. However, note that these "limits" on 
mass and radius assume that the undetected clusters lie on the same 
scaling relations as brighter clusters; if they are in fact relatively 
deficient in X-ray gas this is unlikely to be the case. 

The X-ray image of group 15 shows a clear elongation to the 
north, and is likely a double system. This is a very low surface 
brightness system, from which it is difficult to extract a reliable 
mass profile and hence to estimate 7?5oo . The luminosity and tem- 
perature used here correspond only to the main part of the cluster, 
excluding the northern elongat ion. We have a pplied the Anderson- 
Darling test, as described in H ou et all ([2009^ to all groups, to test 
for non-Gaussianity in the velocity distribution This group is the 
only one that fails the test, with a 99.7 per cent confident detection 
of non-Gaussianity. The spatial distribution of the galaxies is ex- 
tended in the same direction as the X-rays, but this elongation is 
not clearly correlated with the velocity offsets. Thus, in the near- 



infrared data we will not attempt to separate the substructure com- 
ponent. 



2.3 Near-infrared data 

Of the 1 8 clusters in our sample, good-quality near-infrared imag- 
ing was obtained for II, over a period of several years at CTIO, 
AAT and CFHT. For the remaining clusters, the only NIR data 
available is from 2MASS. Table [T] lists the basic properties of our 
clusters and gives the telescope and image size of the NIR observa- 
tions. For the 2MASS clusters, the image size just refers to the area 
used for our analysis — of course the imaging is all-sky. 



2.3.1 CTIO 

Near-infrared data were obtained with the ISPI instrument on 
CTIO, during the nights of May 8-10, 2004. Seeing was ~I" in 
K and two of the nights were fairly cloudy; none were photomet- 
ric. The instrument has a 10.25 x 10.25 arcmin field of view, with 
0.3 arcsec pixels. A 3 x 3 mosaic spanning ^ 30 arcmin was ob- 
tained for eight clusters, but only six were of good enough quality 
to include in the present analysis. Data were reduced in IRAF using 
the CIRRED package, fol lowing the reduction algorithm described 
by lAndreon et alj fcOOOl) . including bias subtraction, flat fielding 
and sky subtraction. The data were calibrated by comparing the 
brightest galaxies with 2MASS, typically involving a few objects 
per field, resulting in a zeropoint accuracy of about ±0.2. The mo- 
saic was combined with the Terapix software SWARP. 



2.3.2 AAT 

In the fall semeste r of 2004 we obtained data with the IRIS 2 
dTinnev et al1l2004l) instrument on the AAT. The field of view of 
this camera is 7.7' x 7.7' with 0.4486" pixels. Six clusters were ob- 
served, by default with a 3 x 3 mosaic, with a resulting coverage of 
23'. We exclude XMM5 from this analysis, for which only the cen- 
tral position was observed. Cluster 1 8 was observed in poor condi- 
tions and several of the pointings were repeated. The final mosaic 
was made with only the best images at each position. 

The data were reduced following the same standard proce- 
dures as for the ISPI data above, but using STARLINK software. The 
astrometric calibratio n was done using the very helpful ASTROM- 
ETRY.NET software tLang et al.ll201ol ; iHogg et alj|2008h . and the 
mosaics were combined with SWARP. Photometric calibration was 
done by comparison of bright galaxies with 2MASS, with similar 
accuracy to our ISPI data. 



2.3.3 CFHT 

Four clusters were observed with WIRCAM on the CFHT during 
the 06A semester. This is a much larger instrument and each cluster 
was observed with a single pointing, covering a 20 arcmin field 
with 0.3 arcsec pixels. Data were reduced by Terapix in Aug 2008. 
The photometric calibration is done by comparison with 2MASS; 
given the larger field the zeropoint is more precise in these fields 
than in our ISPI and IRIS2 data, typically better than ~ 0.1 mag 
rms. 
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Figure 1. The number of galaxies per square degree and per 0.5 magnitude 
bin is siiown for each of the 1 1 clusters with deep NIR observations. The 
black lines and crosses correspond to the CTIO data; red lines and circles 
represent the AAT data; and the blue lines and triangles represent CFHT 
data. 

2.3.4 2MASS 

For the remaining seven clusters in our sample, we do not have deep 
NIR imaging so we use the public 2MASS catalogues. These how- 
ever are much shallower than our data, with a limiting magnitude of 
K < 13.7. It is known that the limiting surface brightness of these 
data is particularly problematic, and the flux from the brightest cD 
galax i es may be signific antly underestimated (e.g. Kochanek et al. 
l200ll;lLin&Mohill2004h . We make no explicit correction for this, 
and simply note that the total luminosities for these clusters may be 
underestimated. 



3 ANALYSIS 

In this section we discuss the data analysis from the eleven clusters 
observed at CTIO, AAT or CFHT. Fluxes were meas ured from the 
reduced data using SEXTRACTOR v2.5.0 ( jSertin & Arnouts 1996 ). 
We use MAG_AUTO as the best estimate of the total magnitude. 
However, star-galaxy separation is done by comparing 4" aper- 
ture magnitudes with the Kron radius. The stellar locus is easily 
identified, separately for each image, and we use this rather than 
CLASS-STAR. This is only relevant for our analysis when cal- 
culating the weights to correct for the spectroscopic sampling rate, 
discussed below. 

We show the number of galaxies per square degree in each 
cluster in Figure [T] The turnover in the counts gives an indication 
of the depth of the observations. The CTIO data are generally shal- 
lowest, with a limit of K ^ 18.0, while the CFHT data are deepest. 

Our NIR photometry is matched to the original 2dF imaging 
catalogues, by searching for the closest match within a 3" aperture. 
The 2dF catalogues are blue-selected, with bj < 19, and this limits 
our sample size. We show a colour-magnitude diagram of our data 
in Figure |2l restricted to galaxies within 1500 km/s of each clus- 
ter. The red sequence of galaxies typical of dense environments is 
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Figure 2. The colour-magnitude relation for observed galaxies with red- 
shifts within 1500 km/s of each cluster. Redshifts from the 2dFGRS were 
obtained only for galaxies with bj < 19; thus our sample becomes incom- 
plete for red galaxies with K > 14.5. 

evident, and the blue selection means our sample becomes incom- 
plete for K > 14.5. This is ~ 3 magnitudes brighter than our K 
detection limit, so there is no additional incompleteness in the NIR 
imaging. 

Brighter than K = 14.5 we calculate a i^— dependent sam- 
pling completeness for each cluster, as shown in Figure [3] This is 
the ratio of the number of galaxies with a redshift to the total num- 
ber detected in the NIR, considering only galaxies within 600kpc of 
each cluster centre. In general, the completeness is close to 80 per 
cent over most of the magnitude range we're interested in, falling 
to 50 per cent al K = 14.5. We weight all galaxies by the inverse 
of this number for all subsequent analysis. In practice the weight is 
calculated separately for each cluster based on a smooth fit to the 
completeness function in the relevant field. 

We use the mean redshift of each cluster to calculate absolute 
magnitudes for all members, using a cosmology of Qm = 0.3, 
JIa ~ 0.7, and h = 0.7. The most distant cluster in our sample is 
al z — 0.096, which corresponds to a luminosity dis tance of 442.0 
Mpc. We correct for Galactic extinction using the 'Schle gel et al.l 
(l998) dust maps (though this is entirely negligib le), and apply a 
k-corr ection of k{z) = — 61og(l + z), following iKochanek et al.l 
( I2OOII) . For our most distant cluster, therefore, we are complete in 
luminosity for Mk < —23.5 . 

First, we present the cumulative luminosity function from the 
combination of all eleven clusters with follow-up NIR data. We in- 
clude all galaxies with redshifts within 1500 km/s of the cluster 
redshift, and within a distance -Rrms from the centr^. We show 
the weighted number of galaxies brighter than a given Alk abso- 
lute magnitude, per cluster in Figure |4] Plotted for comparison are 
Schechter functions with AI^ = —24.3 dBalogh et al. 200 lai) and 
a = —0.5 (solid) or a = —1.0 (dashed). These are not fit to the 
data, but are meant only to guide the eye. The presence of central 

^ For clusters 1, 2, and 5 the maximum radial coverage is defined by the 
extent of the IRIS2 imaging. 
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Figure 3. Top panel: The solid line shows the K magnitude distribution 
of all detected galaxies in our images, while the dotted line shows the dis- 
tribution of galaxies with redshifts from the 2dFGRS. Bottom panel: The 
connected points represent the sampling fraction, which is the fraction of 
ii'-detected galaxies with redshift, as a function of magnitude. Shown here 
is the average selection function; in practice we compute this separately for 
each cluster Only galaxies within 600 kpc of the cluster centre are included. 



Figure 4. The cumulative luminosity function for our eleven clusters with 
deep NIR imaging. This includes only galaxies within 1500km/s of the clus- 
ter redshift, and within Rrms of the centre. The dotted Hne indicates the lu- 
minosity limit con'esponding to Mk > —23.5; at luminosities fainter than 
this our sample is incomplete due to the blue selection of the 2dFGRS. Two 
Schechter functions are shown to guide the eye; they are no t fit to the data. 
Both have = -24.3, taken from iBalogh et alj2001 j) . The sohd line 
has a shallow faint-end slope of a = —0.5, whUe the dashed Une shows 
a = -1.0. 



bright galaxies in most of our clusters leads to an excess relative to 
the Schechte r function, at the brigh t end, as is seen in more massive 
clusters (e.g. IPopesso et ai]l2005bh . Recall that, for Mk > -23.5 
our sample is incomplete, and this limit is indicated by the dotted 
line. 



4 RESULTS 

We now proceed to calculate the total K luminosity of each cluster, 
Lk- We simply sum the luminosity of all galaxies within 1500 km/s 
and i?rms of the cluster centre, brighter than the K = 14.5 limit. 
To correct for galaxies below this limit we model the luminosity 
function as a Schechter function with = — 24.3anda — —1.0. 
As our data reach at least 0.6 mag fainter than Af^ for all clusters, 
this correction is always less than 20%. 

We also include the seven clusters without deeper NIR data 
in our analysis; these are shown as open symbols on the follow- 
ing Figures. For these clusters, the 2MASS limiting magnitude is 
much brighter, K < 13.7, and for three of them the corresponding 
correction for fainter galaxies is larger than a factor of two. 

We do not attempt to measure or correct for intracluster light. 
This remains an important uncertainty in all such work, with some 
claims that a large fraction of the stars in galaxy groups are foun d 
in this component (f Gonzalez et a"Lll2005l ; iMcGee & Baio5ill2010l) . 
For clusters in the mass range of our sample, however, we expect 
the i ntracluster light co ntribution to be less than about 20 per cent 
(e.g. lZibettietal .1120051) . 

The statistical uncertainty on Lk is dominated by the statisti- 
cal uncertainty on Rims, since that quantity determines the radius 




Figure 5. The cumulative i^-band luminosity of our clusters, as a function 
of distance from the centre. This excludes clusters with only 2MASS pho- 
tometry, or for which our NIR coverage does not extend to iJrms- EiTor 
bars show estimated Icr uncertainties arising from Poisson fluctuations in 
the number of galaxies contributing to each bin. The solid line is a fit to the 
data with R/Rrms < 1.0, and has a slope of 0.63. This relation is used to 
calculate how the uncertainty on Lk depends on our uncertainty on -Rrms, 
and to correct the total luminosity of clusters for which our NIR coverage 
does not extend to the radius of interest. 
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Table 3. Derived properties of the stellar population for all clusters in our sample. Throughout the paper, stellar masses are calculated using a universal 
mass-to-light ratio in the K— band, 7 = 0.7. 




Figure 6. Top panel: The dynamical mass of our clusters is shown as a 
function of their total K luminosity, integrated out to Rrms . The horizontal, 
dashed lines indicate the selection limits imposed on the sample, while the 
solid line represents M^y^/Ln = 100. Open circles represent clusters 
with only 2MASS imaging. The green point corresponds to the group 15, an 
outlier on the Lk — Lx relation that will be identified on Figure!?] while the 
red points are undetected in X-rays. Bottom panel: The same, but for A/200 
and the K luminosity within -R200, where these quantities are m easured 
from the velocity dispersion alone, following Ramella et al] j2004h . whose 
data are shown as the crosses without error bars. The solid line represents 
M 20o/L[{ = 100, whi le the dotted line represents Lji oc M*^-^* as found 
bv lRamella et alj j2004l) . 

within which the luminosity is integrated In Figure [5] we show 

Poisson-type uncertainties on the number of cluster members N are not 
really appropriate in this context, as for a specific cluster is a fixed num- 



the cumulative luminosity as a function of r/i?rms, only including 
clusters observed out to -Rrms, with A AT, CTIO or CFHT. The best- 
fit line to the data where R/Rrms < 1 has a slope of 0.63. Thus, 
the statistical uncertainty on Lk is only ~ 0.63A_R/_Rrms which, 
given the typical 10 per cent uncertainty on Rrms, corresponds to a 
~ 6 per cent uncertainty on Lk - In contrast, the typical uncertainty 
on M200 is 20-40 per cent, as it is proportional to (see §|2l. We 
use this relation for Lk{R) to correct the total luminosity of those 
clusters for which NIR coverage only extends out to r < Rrms- 

The total K luminosity of each cluster, integrated to either 
-Rrms, R200 or -R500, is given in Table [3] We also show the total 
number of galaxies with redshifts and NIR data within each radius. 

4.1 Correlation between stellar luminosity and dynamical 
mass 

We have shown above that Lk can easily be measured with a pre- 
cision about five times bet ter than that o f A/dyn, a point that has 
been noted by others (e.g. IPopesso et al]|2005ai) . This makes it a 
very useful indicator of system mass, although of course it is track- 
ing a fundamentally different quantity than Mdyn- In the top panel 
of Figure [6] we show the correlation between these two quantities. 
The clusters were selected to span a factor of only three in A/dyn, 
but they show a factor ~ 10 spread in Lk- Note that the apparent 
lack of correlation is likely a consequence of limited dynamic range 
in A/dyn, together with significant scatter between A/dyn and Lk- 
Most of the clusters are consistent with Mdyn/ Lk = 100, shown 
as the solid line. The green and red points indicate X-ray underlu- 
minous systems, which will be discussed below. Note these most 
of these have Mdyn/ Lk ratios in good agreement with the rest of 

the sample. 

To compare with data from lRamella et ak we compute 

R200 and A/200 in precisely the same way they do, using equa- 
tions [2] and [3] We measure Lk also within R200, and show the 

ber without uncertainty. Statistical uncertainties on the photomefiy, while 
of order 5%, are negligible because they are reduced by \/~N when applied 
to the total luminosity. 
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Figure 7. The bolometric X-ray luminosity of our clusters are compared 
with Lk\ both quantities are computed within Rsop as determined from 
the X-rays. The small, filled circles ms the clusters o f'Lin & Mohj j2004l) . 
with X-ray luminosities taken from lReiprich & Bohr inger 120021), and the 
solid line is a linear fit to these data. The X-ray undetected clusters (red 
points) are represented as upper limits, with Lk measured within iJsoo ^ 
0.4 Mpc. The cluster mai'ked in green is detected in X-ray, but is a clear 
outlier on this plot, with about five times more near-infrared luminosity 
than expected. 



corre lation with M200 in the bottom panel of Figure |6] The data of 
iRame lla et al. (2004) are shown as the small crosses. The solid line 
represents a mass-to-light ratio of M200 / LK{r < R200) = 100, 
while the dotted line shows the relation found by iRamella et al.l 
(|2004) . I/K PC M° Most of the clusters are consistent with the 
data of IRamella et"ai] i2004h . though there are clear outliers. Most 
notable are the two groups (13 and 18) that have stellar masses 
well below that expected from their dynamical mass. The fact that 
they are undetected in X-rays suggests that their dynamical mass 
is significantly overestimated; however, we will show (Fig [7} that 
the lack of X-ray emission is still surprising given their total stellar 
mass. 



4.2 X-ray properties 

In Figure|7]we show the bolometric X-ray luminosity of each clus- 
ter as a function of its total A'— band luminosity. Both quantities are 
computed within R500 as determined from the X-ray emission (Pa- 
per 11). Throughout the paper, both the bolometric luminosity Lx 
and X-ray temperature include contribution from any cool-core. 
This maximizes the scatter and allows us to test for any correla- 
tion between cool-core emission and optical propertie s of the clus- 
ters. W e compare our measurements with the data of iLin & MohJ 
with X-ray luminosities (also uncor rected for cool cores) 
obtained from lReiprich & Bohringei] j2002l) ; the solid line repre- 
sents a least-squares fit to these data. There is considerable scatter 
among our clusters, but most of the X-ray detected clusters are con- 
sistent with those o f iReiprich & BohringeJ ( |2002|) . There is some 
indication that our clusters lie toward the lower edge of the distri- 
bution. This is likely due to the well-known Malmquist bias that 



> 
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Figure 8. The X-ray temperature of our 13 detected clusters are compared 
w ith Lk as measured within -R500 ■ The small, filled circles are the cluster s 
of lLin & Moiiil j2004l) , with X-ray temperatures taken from lHomeil feOOlb . 
The green symbol corresponds to group 15, which is also an outlier in fig- 
ure |2] 



affects X-r ay selected samples (e.g. 'ikebe et al.l l2002l ; IPratt et al] 
l2009l : IVikh linin et al. 2009; Rykoff et al. 2008b). We therefore ex- 
pect our dynamically-selected cluster sample to have a representa- 
tive distribution of Lx- 

One of the systems (group 15, in green) appears to have sig- 
nificantly low Lx for its stellar luminosity. Recall from § 12.21 
that group 15 shows substructure in the X-ray and a significantly 
non-Gaussian velocity distribution. The luminosity and tempera- 
ture used here correspond only to the dominant component. Treat- 
ing the whole system as one results in a temperature increase of 
~ 10 per cent, and a Lx increase of ~ 73%. This does not change 
the status of this system as an outlier on this plot. 

The limits on the five undetected clusters imply that they are 
underluminous in X-ray by up to a factor of 10. Recall that, for 
these clusters the radius i^soo is based on the radius expected for a 
cluster at our iboi detection limit. If the clusters are actually defi- 
cient in X-ray emitting gas, the true value of -R500, and hence Lk, 
could be be larger. This would make the discrepancy worse, since 
the limit on Lx is based on a fixed number of counts in the field of 
view. 

The total X-ray luminosity is known to be a relatively poor 
tracer of mass, but a good tracer of the thermodynamic history 
of the gas (e.g. lBalogh etai]|2006l) . The X-ray temperature, while 
more difficult to measure, is a much better representation of the 
size of the potential. Thus, in Figure|8]we show the correlation be- 
tween temperature and Lk for the 13 clusters in our sample with X- 
ray detections. The temperature measurements are des cribed in Pa- 
er II. Again, these data are compared with data from lLin & Mohr 
2004"), with temperatures (uncorrected for cool cores) from Honiei 
(I2OOI) . Our data generally agree very well with the relation defined 
by the lLin & Mohil i20o4) data, confirming that Lk is a good tracer 
of cluster mass. Interestingly, the detected cluster 15 (in green), that 
is a significant outlier on the Lk — Lx relation, also lies off the 
Lk — T relation in Figure [8] This could indicate that it is a high 
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Figure 9. The correlation between X-ray luminosity and tem perature for 
the 13 detected clusters, compared with the lLin & MohJ |2004) cluster sam- 
ple. The green point indicates group 15, an outlier on the Lfc — Lx rela- 
tion. The fact that it appeal's consistent with the rest of the data shown here 
suggests that the offset is due to unusually high Lfc. 



stellar content, rather than low X-ray luminosity, that makes it un- 
usual. This is also seen in the Lx — T relat ion, shown in F igure[9l 
where our data are compared with those of iHomeil l l200lh . Group 
15 appears consistent with this relation. We again note, however, 
the possib le, small off s et of the whole sample toward lower Lx 
relative to lLin & Moh3 j2004h . indicative of Malmquist bias in X- 
ray flux-limited samples (Paper II). 



4.3 Stellar mass 

Finally, we convert the Tf— luminosity into a stellar mass, assum- 
ing a stellar mass-to-light ratio 7 — 0.7Mq/ Lq jLin & Mohj 
l2004h . In Figure 1101 we compare this with the gas mass within 
i?50o, as measured from the X-ray observations. The lines rep- 
resent different values of the baryon fraction found in stars, / = 
Mstar/ (Mstar + Afgas), of 0.075, 0.12 and 0.2. Our data show 
considerable scatter, but on average our detected clusters are con- 
sistent with ~ 12 per cent of their baryons in stars (neglecting in- 
tracluster light), as typically found for more massive clusters (e.g. 
iBalo gh et al. 2001b). For the five undetected clusters, we show the 
estimated upper limit of A/gas < 1.3 x IO^^Mq, and compute 
the stellar mass within the radius 7?5oo ~ 0.4 Mpc that would be 
consistent with an otherwise normal cluster at our X-ray detection 
limit. These limits imply that > 20 per cent of the expected baryons 
are in the form of stars. If i?5oo is larger than we have estimated, 
for example if their dynamical masses are correct, then Afstar could 
be substantially larger. Thus, our limit on the stellar fraction here is 
robust. These are clearly very different systems from the rest of the 
sample, with little or no associated X-ray emitting gas. 

We next compare the stellar mass with the total mass A/500, 
as estimated from either the X-ray mass profiles (Figure [TTt, or 
the galaxy dynamics (Figure[72l(. In the latter case, we use -R500 = 
V200/500i?200, and M500 = a/200/500 A/200 • The lines in this 
figure represent total stellar fractions of Afstars/A/soo = 0.005, 




Figure 10. The total gas mass within -R500. measured from the X-ray 
emission, is compared with the total stellar mass within the same radius, 
computed from Lx assuming 7 = 0.7Mq / Lq. Lines represent different 
values of / = Mstar/ (A/star + A^gaa) for Comparison. The dotted line 
shows / = 0.12, while the two solid lines show / = 0.075 (top line) and 
/star = 0.2 (bottom). The X-ray undetected clusters are represented as 
red upper limits, with L^ measured within i?5oo ~ 0.4 Mpc. The green 
symbol represents group 15, the outUer on Figure|2] 



0.01 and 0.03. The da ta are compared with lLin & Mohj ( l2004h and 
iGonzalez et alj ( l2007h . For the latter, we exc lude the intracluster 
light component, as neither our data nor thatofLin&Mohr(200' 
includes this. We omit the two clusters in the Gonzalez et al. ( 200' 
sample that are strongly affected by line-of-sight projections; the 
two lowest-mass cl usters remaining in the ir sample have the most 
uncertain R200 and lGonzalez et all ( l2007h exclude them for some 
of their analysis, though we retain them here for the sake of com- 
pleteness. 

In general the agreement with the lLin & Mohj ( |2004|) data is 
good, but our lowest-mass clusters (A^star < 4 x 10^^ A/0) have 
systematically l ower stellar fractions than similar systems from 
iGonzalez et alj ^2007). In particular, none of the X-ray detected 
clusters have total stellar fractions greater than about 3 per cent, 
and most are closer to 1 per ce nt. Our data are in g ood agreement 
with the theoretical prejudice o f'Balogh etal.'teOO^, that the stellar 
fraction in groups and clusters cannot be much larger than observed 
in the most massive systems. 

Some of the five clusters that are undetected in X-rays appear 
to have unusually high stellar fractions; however recall that the ra- 
dius R500 within which the stellar light is integrated is not directly 
measured for these systems. In Figure [72l we use dynamical mea- 
sures of A/500 and i?5oo, which we can measure uniformly for the 
whole sample. Here we see that all our clusters are remarkably uni- 
form, with a total stellar fraction close to ~ 1 per cent in most 
cases. 

Finally, we note that, for the most part, the open symbols on 
Figures[6l410l which indicate the clusters with only 2MASS imag- 
ing, do not appear distinct from the other points on the plots. This 
gives confidence that our results are not significantly biased by the 
lack of deeper imaging in these systems. The obvious exception is 
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Figure 11. The total cluster mass within -Rsoo. as measured from the X- 
ray emission, is compared with the total stellar mass within the same radius, 
computed from Lfc assuming 7 = 0.7Mq/Lq. Lines represent different 
values of / = Mstar/Msoo for comparison. The dotted line shows / = 
0.01, while the two solid lines show / = 0.005 (top line) and /star = 
0.03 (bottom). The X-ray undetected clusters are represented as red points, 
and the green symbol represents group 1 5, the outHer on Figure [T] Small 
cross es are data fror n |Lin & MohJ (2004), and blue triangles are clusters 
from lGonzalez et aP j2007 ), exc luding the intracluster light component, as 
presented in iBalogh et alj 420081) . With the exception of the single outlier, 
all our detected clusters have total stellar fractions less than 3 per cent, 
with an aver age of about 1 per c ent. These stellar fraction s are consistent 
with th ose of lLin & MotiJ J2004 . but lower than found bv lGonzalez et alJ 
i2007h in their lowest-mass systems. 




Figure 12. The same as Figure 1111 but where R500 and A/500 are now 
estimated from the galaxy dynamics, for all clusters in our sample. 
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Figure 13. The velocity distribution for our cluster sample, considering 
galaxies within 1 Mpc of the recomputed centre. The red histogram cor- 
responds to the five undetected systems. The rest of the clusters are in- 
cluded in the black histogram. The smooth lines are Gaussian functions 
with cr = 650km/s, normalized to match the total number of objects in 
each histogram. 

group 15, which is an outlier on most plots. This, however, may 
be related to the fact that it is a clear double-system, and the only 
cluster with a significantly non-Gaussian velocity distribution. 



5 DISCUSSION 

Of the 18 clusters observed, five are undetected in X-ray, and one 
(group 15) appears to have an anomalously high stellar content 
given its X-ray properties. The five non-detections have a stellar 
mass which is ~ 1 per cent or less of their total dynamical mass, 
consistent with most other clusters. However, the strong limits on 
their X-ray luminosity imply that there is little or no associated X- 
ray emitting gas. For these systems, either their virialised mass is 
M < Mdyn, or they are extremely underluminous due to excess 
heating or cooling. 

Of these, the least likely possibility is that all the gas has 
cooled to form stars. Figure [T2I suggests that, if anything, the ra- 
tio of stellar mass to dynamical mass in these clusters is lower than 
that of the normal systems. Normal clusters have ~ 7 times more 
gas mass than stellar mass (Figure fTOt; if all that gas formed stars 
we would expect to see stellar luminosities ~ 7 times higher than 
normal clusters of the same dynamical mass, and this is clearly not 
the case. 

On the other hand it is difficult to rule out the possibility 
that the gas in the underluminous clusters has been heated or 
expelled, as predicted by some models (e.g. lBoweretal.ll2008l : 
[McCarthy et al.l2010j) . In this case of course there should be some 
gas present, but at low density and low surface brightness. The only 
way to definitively distinguish between this and the following sce- 
nario is through X-ray observations deep enough to detect either 
this hot, diffuse gas or the gas associated with individual galaxies. 

This leaves the possibility that the undetected clusters are 
unrelaxed, and have not yet reached virial equilibrium (e.g. 
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Figure 14. For each cluster, we show the cluster X-ray luminosity as a 
function of the ii'-band luminosity of the brightest cluster member (BCG). 
Point symbols are as in other plots, with the green symbol for the outlier 
group 15, and red symbols for undetected clusters. Clusters are numbered 
according to their entry in Table [T] Open symbols represent clusters for 
which the near-infrared imaging comes from 2MASS. The lack of correla- 
tion here demonstrates that the clusters which are underluminous in X-ray 
host BCGs that are just as luminous as the rest of the sample. 



IPopesso et al.l2007l ; lDietrich et alj26o9l : lRasmussen et alj2006h . It 

is even possible that some of these systems are chance projections 
, in which case the appropriate "virialised mass" to consider is just 
that of the dominant galaxy. To shed light on this issue we now 
investigate this using the spatial and dynamic information avail- 
able for each cluster. Since these results are somewhat sensitive to 
the choice of centre and radius, we recompute these quantities for 
all the following analysis, including galaxies with redshifts from 
surveys other than the original 2dFGRS. Specifically, we select all 
galaxies with redshifts, within 1.5 Mpc of the original centre and 
1500 km/s of the original cluster redshift. We adopt the central po- 
sition and redshift as the geometric mean of these quantities, un- 
weighted for selection or luminosity. 

In Figure [T3] we show the velocity distributions for the "nor- 
mal" clusters (black line), including only galaxies within 1 Mpc of 
the recomputed centre. We compare this with the five undetected 
clusters, as the red line. Overplotted for comparison are Gaussian 
functions with a velocity dispersion of 650km/s, for comparison. 
A Kolmogorov-Smimov test shows that the velocity distributions 
of the undetected clusters are consistent with that of the rest of the 
sample, with a 13 per cent probability of either being drawn from 
the same parent distribution. Thus, given the limited statistics, we 
find no evidence for differences in the dynamical state of the under- 
luminous systems. As described in § 12.21 we h ave also compute d 
the Anderson-Darling statistic, as described by IHou et al] j2009l) , 
for each group. Only group 15 shows significant non-Gaussianity 
in the velocity distribution, at > 95 per cent confidence. 

Other indications of a relaxed cluster co uld be the presence o f 
a dominant galaxy near the X-ray centre (e.g. iDariush et al1l2010h . 
and an approximately spherical, centrally-concentrated galaxy dis- 
tribution. In many of the clusters, there is an obvious dominant 
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Figure 15. For each cluster, we show the luminosity ratio of the bright- 
est and third-brightest galaxies (in K), as a function of the distance of the 
brightest galaxy from the recomputed cluster centre. Symbols are as in Fig- 
ure [14] 



galaxy, large and luminous, near the centre of the X-ray image. To 
quantify this, we have selected all cluster members within 700 km/s 
of the mean redshift and 500 kpc of the recomputed centre, and 
identified the most luminous (in K) as the BCG. Figure [14] shows 
the luminosity of each of these BCGs, as a function of the host 
cluster Lx- There is little correlation here, although we note that 
three of the X-ray undetected clusters have BCGs that are among 
the least luminous in the sample (unfortunately, all three have only 
shallow 2MASS imaging, and thus the total luminosities may be 
significantly underestimated). It is also interesting that group 15, 
which has a high stellar mass given its Lx, has one of the most 
massive BCGs in the sample. 

Next we calculated the luminosity ratio between the first- and 
third-ranked galaxy, Lk.is- In Figure [TS] we show this as a func- 
tion of A_R, which is the distance in Mpc between the brightest 
cluster galaxy (BCG) and the geometric centre of the cluster, re- 
computed as above. Interestingly, the distribution of the X-ray un- 
detected clusters (and also the outlier group 15) is distinct from 
most of the "normal" systems, in the sense that their BCG is at least 
250 kpc from the centre. However, we note that, of the 12 otherwise 
"normal" systems, only about half have a dominant (Lk ,13 > 3), 
centrally located (Ai? < 250 kpc) BCG. 

We now attempt to quantify the spatial distribution of the most 
luminous galaxies, i.e. those that are at most 0.6 mag fainter than 
M^, so that we are equally deep in all clusters. The undetected 
cluster 18 only has three galaxies above this limit, so we omit it 
from the following analysis. We calculate the concentration as the 
fraction of such cluster members within 0.5i?rms. For the elonga- 
tion, we first perform a least-squares regression analysis to find the 
principle axis of each cluster on the sky; then we calculate the rms 
dispersion perpendicular to and parallel to this axis. The elonga- 
tion is the ratio of the two values, always defined as the larger 
divided by the smaller so the ratio is greater than unity. The re- 
sults are shown in Figure [161 vvhere the points are colour-coded 
as before. Error bars are computed using a jackknife resampling. 
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Figure 16. For each cluster (except cluster 18, which has too few members 
for this analysis), we show its elongation as a function of concentration. 
Both quantities are computed for a subset of cluster members, brighter than 
+ 0.6 and with velocities within twice the cluster velocity dispersion. 
The concentration is the fraction of galaxies that lie within 0.5i?rms, while 
the elongation is the ratio of the nns values about the directions perpendic- 
ular to and parallel to the least-squares fit congelation in spatial coordinates. 
Symbols are the same as figure 1151 Three of the undetected clusters show 
low values of concentration, and two of these are also highly elongated. 



Interestingly, most of the X-ray underluminous systems again ap- 
pear separated from the majority of the population, as either low- 
concentration or highly elongated clusters. Only cluster 17 lies in 
a region of the plane occupied by the majority of normal clusters. 
Note that cluster 16 is highly elongated; from Figure[T5]we see that 
it has a distinctly dominated galaxy, but located 350 kpc from the 
centre. This may indicate a merging or otherwise unrelaxed sys- 
tem. Cluster 18, which is the only other X-ray undetected system 
with a dominant, centrally-located galaxy, has too few members to 
measure either quantity shown here with adequate precision. Again, 
however, there are examples of clusters (3 and 9) with normal X-ray 
properties and equally low concentrations. 

This result needs to be approached with some caution, as there 
are multiple parameters at work here, related to the magnitude and 
velocity selections, the choice of centre, and the definition of con- 
centration. Nonetheless, we tentatively conclude that it seems likely 
dynamical age plays some role in the X-ray luminosity of a given 
cluster. All the clusters that appear relaxed in the optical - with a 
dominant central galaxy (within 250kpc of the cluster centre, and 
at least three times brighter than the third-ranked galaxy), a cen- 
trally concentrated (> 0.25) galaxy population, and a spatial axis 
ratio of less than two - show normal X-ray propertied There is no 
evidence for excess heating or cooling in these systems. However, 
this does not mean that X-ray selection returns a sample of relaxed 
clusters; of the twelve clusters with normal X-ray properties (i.e. 
excluding the outlier 15 and the undetected clusters), seven show 



^ The possible exception is group 18, which has a dominant, near-central 
galaxy but for which we are unable to measure elongation and concentra- 
tion. 



signs of dynamical youth. In other words, neither optical- nor X-ray 
selection returns a sample of dynamically relaxed clusters, which 
is certainly not a surprise. It is worth emphasizing though that any 
signs of dynamical disturbance are fairly subtle. With the possible 
exception of group 15, none of our clusters look like strongly merg- 
ing systems either in the X-ray or optical images, and the velocity 
distributions (based on few members) do not show obvious signs of 
nonvirialisation. 

Thus, with the present data, we can plausibly argue that the 
undetected systems might be unvirialised and that this is the ori- 
gin of the low X-ray emission. On the other hand, the differences 
between the luminous and under-luminous systems are remarkably 
subtle and we cannot rule out that the difference is due to a sub- 
stantial amount of ejected gas, as predict ed by some models (e.g. 
iBower et al.ll2008l;lMcCarthv et al.ll2(}T0al) . 



6 CONCLUSIONS 



We hav e presented a sample of 18 clusters selected from lEke et al.l 
l l2004a l) based on their dynamical mass alone, and followed up six- 
teen of them with pointed Chandra and XMM observations. The 
clusters were selected based on optical properties alone, to lie in a 



narrow range of dynamical mass 3 x 10 
We summarize our findings as follows: 



< M/Mq < 6 X 10" 



• Of the 18 clusters studied, five are undetected in X-ray emis- 
sion. The limits on their X-ray luminosity are significantly below 
the luminosity that would be expected given their stellar mass. The 
rest lie within the scatter defined by X-ray selected samples, but to- 
ward the low-Lx end of the scatt er. This is in go od agreement wit h 
previous studies b y e.g. Bower et al. (1997), Rv koff et al.l ( l2008ah . 
iHicks et all ( |2008|) and lRasmussen et all a2006i) . 

• Stars make up less than 3 per cent of the total mass in all of 
our clusters, with an average of about 1 per cent. The fraction of 
detectable baryons in the form of stars is about 12 per cent on av- 
erage, with a range of approximately 7—20 per cent, though the 
clusters that are undetected in X-ray emission have implied frac- 
tions of > 20 per cent. 

• The undetected clusters have velocity distributions that are not 
significantly different from the rest of the sample. However, they all 
either lack a central, dominant galaxy, or show spatial distributions 
that are of low concentration or high elongation, relative to most 
of the "norma l " clus ters. Similar conclu s ions h ave been reached by 
IPopesso et all ilOOH) and loietrich et al.l ( |2009|) . 

We conclude that the redshift-selection of low-mass clusters 
returns a heterogeneous sample. Most of those systems that are de- 
tected in X-ray have gas properties that are consistent with X- 
ray selected samples, though somewhat underluminous on aver- 
age. The main difference is that a fraction (6/18 in our case) have 
surprisingly low X-ray luminosity, given their stellar content. It 
remains an open question, whether this is because such clusters 
are dynamically young (unvirialised) or because the hot gas has 
been expelled or rarefied due to energetic processes associated with 
galaxy formation. 
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